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Tuning the trip to KPZ asymptopia 

The paradigmatic stochastic partial differential equa- 
tion proposed by Kardar, Parisi and Zhang (KPZ) has, 
in the past decade, become the cornerstone in our under- 
standing of kinetic roughening phenomena; mathemati- 
cal links to directed polymers in random media (DPRM) , 
a baby version version of the spin-glass problem exactly 
soluble in d = 1, as well as driven diffusive systems, do- 
ing much to intensify community interest in the strong- 
coupling, nonperturbative scaling behavior of this fun- 
damental equation Q. Central to the KPZ dogma is 
the implicit assumption of universality; i.e., the equa- 
tion provides a bare-bones, but catch-all continuum de- 
scription that transcends all irrelevant microscopic de- 
tails, particularly those associated with the lattice. In- 
deed, within the immediate context of noncquilibrium 
stochastic growth, the KPZ equation reigns supreme- its 
very broad universality has been numerically established 
well beyond characteristic scaling exponents, to include 
finite-size corrections, amplitude ratios and full proba- 
bility distributions This Comment arises from a re- 
cent dialog concerning the suggestion of possible KPZ 
nonuniversality, initiated by Newman and Swift [||, who 
numerically integrate the KPZ equation in the strong- 
coupling limit, where their prescription reduces to a zero- 
temperature DPRM analysis. For gaussian noise dis- 
tributions, they readily retrieve accepted values for the 
KPZ energy fluctuation exponent f3 in d = 1 — 4. For 
bounded, symmetric, and innocuous distributions of the 
form: P(e) ~ (1 — |e|) a , however, they notice increasingly 
greater deviations as a — > — 1. In fact, the discrepancies 
grow more severe in higher dimensionalities. This be- 
havior is easily seen as a consequence of the close blood 
relation between the DPRM and its near cousin, directed 
percolation (DP), as emphasized by several authors |||| , 
who point out the important role of falling percolation 
thresholds with increasing dimensionality. Our purpose 
here is clarify the essential features of this tale, insur- 
ing that crucial bits of evidence remain atop the table. 
In this manner, one might hope to discover the optimal 
road to KPZ asymptopia. 

Much can be gleaned, already, using truncated gaus- 
sian noise, restricting ourselves to either positive or neg- 
ative bond energies, thereby breaking the symmetry in- 
trinsic to the full distribution. Interestingly, these demi- 
gaussians yield strikingly different scaling behaviors in 
the zero-temperature DPRM global minimization proce- 
dure. For a positive-bond demi-gaussian, P(e > 0) ~ 
e~ e I 2 , we obtain an anomalous exponent that is much 
too small, though entirely reminiscent of Newman-Swift 
results; e.g., in a modest, but highly representative, sim- 
ulation of the 2+1 dimensional DPRM, involving 10 5 re- 
alizations of the random energy landscape and directed 
paths of length 10 2 steps, the extracted fluctuation ex- 
ponent, (3 w 0.116, is well below the accepted value 
0.24+ ||. Naively, this result appears counterintuitive, 
since a positive demi-gaussian is strongly biased low, pos- 
sessing a relatively high percentage of favorable, low en- 
ergy bonds attractive to the directed polymer. Yet, it is 
precisely this abundance of riches which is responsible for 



a near ground-state degeneracy that greatly delays the 
DPRM's true approach to asymptopia. By contrast, for a 
negative-bond demi-gaussian, P(e < 0) ~ e~ £ "/ 2 , which 
is biased high, we retrieve excellent KPZ scaling; i.e., a 
simulation, otherwise identical to that just mentioned, 
now produces an exponent « 0.245, even for these rather 
short optimal paths! The relative scarcity of particu- 
larly favorable bonds permits, somewhat unexpectedly, 
an accelerated convergence to asymptopia. Surprisingly, 
a simulation employing the full gaussian actually fares 
less well, producing an exponent rs 0.216 under the same 
conditions. 

Splitting the Newman-Swift bimodal (a = —1) distri- 
bution into two pieces unearths, of course, the very same 
pathology- an immediate DPRM vestigial link to directed 
percolation (||| • Biasing high, with positive bond ener- 
gies only, generates excellent KPZ scaling; as impressive, 
in fact, as the previous demi-gaussian! By contrast, bi- 
asing low yields a disasterously small effective exponent. 
The dangers inherent in the Newman-Swift class of dis- 
tributions, although most apparent for a < 0, are also 
clearly present for a > 0, the essential difference being 
one of degree, rather than kind, since all such bare prob- 
ability distributions lie within the DPRM basin of at- 
traction, albeit at different distances from the inevitable 
fixed point function. Indeed, the kernel of the Newman- 
Swift numerical mystery well predates their paper, and 
can be found originally in the work of Kim, Bray and 
Moore , wherein it is revealed that a flat distribution 
(a = 0), while the gaussian's peer within the 1+1 DPRM 
context, seriously underperforms the latter in 2+1 di- 
mensions, for which the DP threshold has fallen below 
one-half fe|. Furthermore, earlier work of Derrida and 
Griffiths H had indicated that a sufficiently asymmet- 
ric, discrete bare probability distribution would land the 
model well outside the KPZ universality class, controlled 
by the DP fixed point. Problems inevitably arise when- 
ever there is too much integrated weight in the low end 
of the chosen DPRM bond distribution 0. Such is the 
price paid for a surprising, but deep blood connection to 
directed percolation. 
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2+1 Wedge DPRM 

nt-100, nr= 10**5, full & demi-gaussians 
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■ demi- gaussian, biased high 
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